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POROUS MEDIUM FLOW WITH BOTH A FRACTIONAL 
POTENTIAL PRESSURE AND FRACTIONAL TIME 
DERIVATIVE 

MARK ALLEN, LUIS CAFFARELLI, AND ALEXIS VASSEUR 


Abstract. We study a porous medium equation with right hand side. The 
operator has nonlocal diffusion effects given by an inverse fractional Laplacian 
operator. The derivative in time is also fractional of Caputo-type and which 
takes into account “memory”. The precise model is 

T>fU — div(-!i(—A)~‘^u) = /, 0 < (T < 1/2. 

We pose the problem over {t S M”*", x G R"} with nonnegative initial data 
^(0, i) > 0 as well as right hand side / > 0. We first prove existence for weak 
solutions when /, tt(0, x) have exponential decay at infinity. Our main result 
is Holder continuity for such weak solutions. 


1. Introduction 

In this paper we study both existence and regularity for solutions to a porous 
medium equation. The pressure is related to the density via a nonlocal operator. 
This diffusion takes into account long-range effects. The time derivative is nonlocal 
and fractional and therefore takes into account the past. In the typical derivation 
of the porous medium equation (see [M]) the equation one considers is 

dtu + div(vu) = 0, 

with u{t, x) > 0. By Darcy’s law in a porous medium v = — Vp arises as a potential 
where p is the pressure. According to a state law p = f{u). In our case we consider 
a potential which takes into account long range interactions, namely p = {—A)~'^u. 
A porous medium equation with a pressure of this type 

(1.1) dtU = div{u{—A)~'^u) 

has been recently studied. For 0 < cr < 1 with a ^ 1/2, existence of solutions was 
shown in [S] while regularity and further existence properites were studied in [3]. 
Uniqueness for the range 1/2 < tr < 1 was shown in [T7]. Another model of the 
porous medium equation 

— diy{K{u)Du) = f 

was introduced by Caputo in [^ . In the above equation Df is the Caputo derivative 
and the diffusion is local. Solvability for a more general equation was recently 
studied in m- The fractional derivative takes into account models in which there 
is “memory”. The Caputo derivative has also been recently shown (see [7], [8] ) to 
be effective in modeling problems in plasma transport. See also [12] and [16] for 
further models that utilize fractional equations in both space and time to account 
for long-range interactions as well as the past. 
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The specific equation we study is 
(1.2) V^uit, x) — div (uV(—A)“'^u) = f(t, x). 

The operator is of Caputo-type and is defined by 

a 


V?u := 


[u(t, x) — u{s^ x)]K{t^ s, x) ds. 


r(l-a) J_, 

When K{t, s, x) = {t — this is exactly the Caputo derivative - see Section[2] 

- which we denote by Df. We assume the following bounds on the kernel K 
1 . Of ^ A 


(1.3) 


< 


-K{t, s, x) < 


A(t-s)i+“ - r(l-a) " ’ '-(t-s)i+“ 

Our kernel in time then can be thought of as having “bounded, measureable coef¬ 
ficients” . We also require the following relation on the kernel 


(1.4) K{t,t — s) = K{t + s,t). 

The relation (d allows us to give a weak - in space and in time - formulation of 
d- This weak formulation is given in Section [2l 

In this paper we also restrict ourselves to the range 0<cr<l/2. In[4] use of a 
transport term was made to work in the range 1/2 < cr < 1. We have not yet found 
the correct manner in which to prove our results for 1 /2 < ct < 1 when dealing with 
the nonlocal fractional time derivative 2 ?/. 


1.1. Accounting for the Past. Nonlocal equations are effective in taking into 
account long-range interactions and taking into account the past. However, the 
nonlocal aspect of the equation provides both advantages and disadvantages in 
studying local aspects of the equation. One advantage is that there is a relation 
between two points built into the equation. Indeed, we utilize two nonlocal terms 
that are not present in the classical porous medium equation to prove Lemma 
15.31 One disadvantage of nonlocal equations is that when rescaling of the form 
v{t,x) = Au{Bt,Cx), the far away portions of u cannot be discarded, and hence 
V begins to build up a “tail”. Consequently, the usual test function {u — fc)+ or 
F((u — fc)+) for some function F and a constant k is often insufficient. One must 
instead consider F((tt—</>)+) where is constant close by but has some “tail” growth 
at infinity. This difficulty of course presents itself with the Caputo derivative. One 
issue becomes immediately apparent. If we choose F{{u — </)+) as a test function, 
then 






/')+)A“((« 


</)+ 



0)+)24“</ dt. 


{u — (j))-) dt 


The second term will no longer be identically zero if 4> is not constant. When 
using energy methods, this second term can be treated as part of the right hand 
side, and hence it becomes natural to consider an equation of the form (O with 
a right hand side. The main challenge with accomodating a nonzero right hand 
side is that the natural test function Inu used in [5] and [4] is no longer available 
since the function u can evaluate zero. Indeed if the initial data for a solution is 
compactly supported, then the solution is compactly supported for every time t > 0 
(see Remark EH). We choose as our basic test function u"' for 7 > 0. For a small 
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we will have to choose 7 small. We then can accomodate a right hand side as well 
as avoiding delicate integrability issues involved when using In u as a test function. 
Using careful analysis, it is still most likely possible to utilize In u as a test function 
for our equation CJ) with zero right hand side, but we find it more convenient to 
use and prove the stronger result that includes a right hand side. Our method 
using should also work for the equation HU to be able to prove existence and 
regularity with a right hand side. One benefit of accomodating a right hand side 
in L°° is that we obtain immediately regularity up to the initial time for smooth 
initial data, see Theorem ll.41 

1.2. Overview of the Main results. We will prove our results for a class of weak 
solutions (12.4|) later formulated in Section [2] Our first main result is existence. We 
use an approximating scheme as in [5] as well as discretizing in time as in [1]. We 
prove 

Theorem 1.1. Let 0 < uq{x), f{t,x) < for some A > 0. Then there exists 

a solution u to (in in ( 0 ,oo) that has initial data u{0,x) = uo(x). 

Remark 1.2. Our constructions are made via recursion over a finite time interval 
(OjTi). Since our constructions are made via recursion, if T 2 = mTi for m G N, 
it is immediate that if Ui is the solution constructed on (0,Ti), then U 2 = ui on 

(o,ri). 

Remark 1.3. For technical reasons seen in the proof of Lemma 14.41 when n = 1 we 
make the further restriction 0 < cr < 1/4. 

The main result of the paper is an interior Holder regularity result. As expected 
the Holder norm will depend on the distance from the interior domain to the initial 
time to- However, if we assume the intial data uq is regular enough - say for intance 
- then we obtain regularity up to the initial time. This is a benefit of allowing 
a right hand side. By extending the values of our solution u{t, x) = u(0, x) for 
t < 0, we satisfy (12.41) on (— 00 , 00 ) x K” with a right hand side in L°°. The right 
hand side / for t < 0 will not necessarily satisfy / > 0 ; however, this nonnegativity 
assumption on / was only necessary to guarantee the existence of a solution u > 0 . 
It is not a necessary assumption to prove regularity. From Remark ll.2l the solution 
constructed on (— 00 , T) will agree with the original solution over the interval (0, T). 

Theorem 1.4. Let u be a solution to (1^ obtained via approximation from Theo- 
rem \l.l\ on [0.r]xR" withO < un(x). f(t.x) < Assume also uq € C'^. Then 

u is continuous on [0,r] x R"- for some exponent (3 depending on a,A,n,a - 
with a constant that depends on the L°° norm of u and f and norm of uq . 

1.3. Future Directions. We prove existence and regularity for solutions obtained 
via limiting approximations. In this paper we do not address the issue of uniqueness. 
As mentioned earlier, uniqueness for (HU for the range 1/2 < cr < 1 was shown in 
m- The issue of uniqueness for HU is not trivial because of the nonlinear aspect 
of the equation as well as the lack of a comparison principle. The equation (12.41) 
which we consider also should present new difficulties because of the weak/very- 
weak formulation in time as well as the minimal “bounded, measurable” assumption 
HU on the kernel K{t, s, x). An interesting problem would be to then address the 
issue of uniqueness for solutions of dMl). 

Theorems II.II and II.41 can most likely be further refined by making less assump¬ 
tions on u( 0 ,a;), assuming a right hand side f G L^ as was done for a similar 
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problem in |15] . and proving the estimates uniform as cr —^ 0 and recoving Holder 
continuity for the local diffusion problem. Also, as mentioned earlier the Theorems 
can be improved to include the range 1/2 < cr < 1 . 

Finally, just like in the local porous medium equation |14j as well as for (HD, 
the equation (|2.4I) has the property of finit propagation, see Remark l4.6l Therefore, 
it is of interest to study the free boundary d{u{t,x) > 0 }. 

1.4. Outline. The outline of this paper will be as follows. In Section [5] we state 
basic results for the Caputo derivative. We also give the weak formulation of the 
equation we study. In Section [3] we state some results for the discretized version of 
Vf that we will use to prove the existence of solutions. In Section 2] we follow the 
method approximation method and use the estimates from [5] combined with the 
method of discretization and the estimates presented in [1] to prove existence. In 
Section [5] we state the main Lemmas that we will need to be able to prove Holder 
regularity. In Section [ 6 ] we prove the most technically difficult Lemma 15.11 of the 
paper. This Lemma l5.1l most directly handles the degenerate nature of the problem. 
In Section [3 we prove an analogue of l5.ll In Section [ 8 ] we prove the final Lemmas 
we need which give a one-sided decrease in oscillation from above. The one-sided 
decrease in oscillation combined with Lemma IQ is enough to prove the Holder 
regularity and this is explained in Section [31 

1.5. Notation. We list here the notation that will be used consistently throughout 
the paper. The following letters are fixed throughout the paper and always refer 
to: 

• a - the order of the Caputo derivative. 

• a - the order of inverse fractional Laplacian (—A)“°’. We use a for the 
order because s will always be a variable for time. 

• a - the initial time for which our equation is defined. 

• Uf - the Caputo derivative as defined in Section [3] 

• 2?/ - the Caputo-type fractional derivative with “bounded, measurable” 
coefficients with bounds dESD and relation (El- 

• A the constant appearing in (11.311 . 

• - the discretized version of 2?/ as defined in (13.111 

• e - will always refer to the time length of the discrete approximations as 
defined in Section [3] 

• n - will always refer to the space dimension. 

• Tm - the parabolic cylinder (—m, 0) x Bm- 

• - the fractional Sobolev space as defined in [9]. 

• u± - the positive and negative parts respectively so that u = u+ — U-. 

• u - the extension u{t) = u{ej) for ej — 1 < t, ej. 


2. Caputo Derivative 


In this section we state various properties of the Caputo derivative that will be 
useful. The Caputo derivative for 0 < a < 1 is defined by 
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By using integration by parts we have 


( 2 . 1 ) 


r(l - a) aD^t) 


u{t) — u{a) A* u{t) — u{s) 

(t-a)“ (i - s)“+i 


ds. 


For the remainder of the paper we will drop the subscript a when the initial point 
is understood. We now recall some properties of the Caputo derivative that were 
proven in [T]. 

For a function g(t) defined on it is advantageous to define g{t) for t < a. 

Then we have the formulation 


aD^g{t) 


ooD^g{t) 


« /■* 9(.t)-9is) 

r(l - o) J-oo {t - s)i+“ 


utilized in [T]. (See also [5] for properties of this one-sided nonlocal derivative.) This 
looks very similar to (—A)“ except the integration only occurs for s < t. In this 
manner the Caputo derivative retains directional derivative behavior while at the 
same time sharing certain properties with (—A)“. This is perhaps best illustrated 
by the following integration by parts formula for the Caputo derivative 


Proposition 2.1. Let g,h G C^(a,T). Then 


gDth + hD^g = f g{t)h{t) 

J a 


1 1 
-f 


( 2 . 2 ) 


a 


_(T-t)°‘ (t-a)“_ 

r [git)-g{s)][h{t)-h{s)] 

la Ja {t-sy+- 

g{t)h[a) + h{t)g{<^ 


dt 

ds dt 


Ja it - o)“ 

Formula (EH) is based on the following formal computation 

r 9is) 


la Ja (t-s)l+“ 1 ^ ,, ,, 


T pt 


J I 


{t - s)i+“ 

9it) 

{t - s)i+“ 


1 +a 


ds dt — 


pT rt 


9 it) 


pT pT 


'a J s 
pT pT 


-[ .(/ 


[t - s)i+“ 
ds 


ds dt — 


J a J t 
pT-t 


it-s) 

gjs) 

it - s)i+“ 

gjt) 

is-t) 


l+ct 


ds dt 

dt ds 

ds dt 


ds 


oH-a 


ol + Ct 


dt 


In the above computation to utilize the cancellation we only need a kernel K(t, s) 
satsifying 

K{t,t — s) = K{t -\- s,t). 

To make the above computation rigorous we will use the discretization in Section 
|3l An alternative, equivalent integration by parts formula is to extend git) = g(a) 
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for t < a. Then for any h € with h{t) = 0 for any f < 6 for some b we have 


h{t)D°g{t) dt = Co 


(2.3) 


r2t-T 


' —OO — oo 
pT rt 


[gjt) - 9 is)][Ht) - h{s)] 
{t - s)i+“ 

g{t)h(t) 


ds dt 


/-OO [t - 


ds dt 


g{t)Dfh{t) dt, 


with Cq = ar(l — a) ^. Now (12.21) and (12.31) will both imply each other, so (12.31) 
is an alternative way of handling the initial condition g{a). Furthermore, both 
(1^ and (|2.3I) are weak formulations for the Caputo derivative that only require 
that g G T)). (12.31) will work for any kernel IC(t, s) satisfying the relation 

(HH). In view of (12.3|) we now give the exact formulation of our weak solutions. We 
assume the bounds (O and the relation (HU on the kernel K{t, s, x). For smooth 
initial data uq G C^, we assign u{t, x) = u{a, x) for t < a. Then as stated earlier in 
the Introduction, for t < a, a solution u will have right hand side 


div(uo(—A) '^uo) G L°°. 

We say that u is a weak solution if for any (jj G C^{—oo, T) x K" we have 


/R" J - 


[u{t, x) — u{s, x)][4>{t, x) — 4>{s, x)]K{t, s, x) ds dt dx 


OO J —OO 


r2t-T 


+ / / / u{t,x)4>{t,x)K(t, s,x) ds dt dx 

JmV- j — OO J — OO 

(2.4) — [ [ u{t,x)'D^cj){t,x) dt dx 

j — OO 

+ f f V(j){t,x)u{t,x)'V{—A)~’^udxdt 
J — OO 

= [ [ f{t,x)(t){t,x). 

J —OO ^ Br 

We will also utilize a fractional Sobolev norm that arises from the fractional 
derivative. 


Lemma 2.2. Let u be defined on [a, T], We have for two constants ci,C 2 depending 
on a, \T — a\ 


||u||^.^(a,T) < ci||u||^„/2(a,T) 


< C 2 



|M(t) -u(s)P 
\t - s|i+“ 


ds dt + 



(T-a)-) 


The following estimate will be needed for our choice of cut-off functions. 


Lemma 2.3. Let 

h{f) := max{|t|‘" — 1,0} 


with v < a. Then 

Vfh > -Cu,a,K 

for t G R. Here, is a constant depending only on a, u, A. 
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Finally, we point out that ii g = g+ — g- the positive and negative parts respec¬ 
tively, then 


(2.5) 


[ g±{t)V^g^{t)>0. 
J a 


3. Discretization in time 


To prove existence of solutions to (12.41) we will discretize in time. The discretiza¬ 
tion also allows us to make the computations involving the fractional derivative 
rigorous. This section contains properties of a discrete fractional derivative which 
we will utilize. 

For future reference we denote the discrete Fractional derivative with kernel K 
as 

(3.1) V‘^u{a + ej) ■.= e ^ [u{a + ej) — u{a + ei)\K{ej, ei) 

The following is the discretized argument for the cancellation that appears in 
the formal computation of the version of (EB that we will need. 


Lemma 3.1. Assume g{a) = 0 and define g(t) := 0 for t < a. Assume relation 
(11.41) on K. Then 

e‘^^Yl[sia + ej)-g{a + ei)]K{a + ej,a + ei) = e‘^Y^ X! 

j^k i<.j j^k i<.2j — k—l 


Let g{f) := g{ej) for ej — 1 < t < ej. If g >0, then there exists c depending only 
on a, A such that if e <1, then 


+ ej) - 9{a + ei)]K{a + ej, a + ei) 

j<k i<j 



(T - t)“ 


dt. 


Proof. For notational simplicity we assume a = 0. 

e^X15Z[5(ej) -5(eO]^(ej.eO 

j<ki<j 

g{ej)K{ej, ei) - ^^EE g{ei)K{ej, ei) 

j<k i<j j<k i<j 

-^EE g{ej)K{ej, ei) - ^^E E g{ei)K{ej, ei) 

j<k i<j i<.k i<j<k 

-^EE g{ej)K{ej, ei) - ^^E E g{ej)K{ei,ej) 

j<k i<j j<.k j<.i<k 

= ^^E E 9{(^j)K{(^j,ei) 

j<k i<.2j — k — l 

The second inequality follows from the estimates in [1]. □ 


Lemma o combined with the estimates in [1] can be used to show. 

Lemma 3.2. Let u(0) = 0 and assume u > 0. For fixed 0 < e < 1, let u be the 
extension defined as in Lemma \3.1\ Let K satisfy conditions (11.31) and (11.41) . Then 
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there exists two constants ci, C 2 depending only on a such that 

) 1 —a 

= 

where c depends on a and A. 


u{ej)[u{€j)-u{ei)]K{€j,ei) > c||u ||^,,/2 > c||m||^,,/ 2 C | f 

j<k i<j \ “ 


This next lemma is analogous to Lemma 12.31 and was shown in [T] . 

Lemma 3.3. Let h be as in Lemma \2.3\ . Then for 0 < e < 1 there exists c^^a 
depending on a and v but independent of a such that 

'T>fh{t) > -c^^a 


for t S eZ and a < t < 0. 

This last estimate we will use often 


Lemma 3.4. Let F be a convex function with F" >'y,F'> 0,F(0) = 0. Assume 
g ^ 0,g{a) = 0. Then there exists c depending on a, A such that 


j<k 


F{g{ej)) 


7 2 

C —€ 
2 


EE 

j<k e<j 


(e(j-*))!+“ 


Proof. Since F is convex, 

-P'(5(ej))b(ej) - g(e0] > F{g{ej)) - F{g{ei)) 
The result then follows from applying Lemma l3.II 




□ 


Finally we point out that if g is a limit of ge which are discretized problems with 
the assumptions in Lemma 13.41 it follows that 


(3.2) 


F{ 9 {t))V^g{t) dt>c 


r Fjgjt)) 

L (T-a)“ 

4. Existence 


c- 


[9{t)-g{s)r^ 

{t - s)i+“ 


ds dt. 


In this section we prove the existence of weak solutions following the construction 
given in [^. We will also discretize in time. We first consider a smooth approxima¬ 
tion of the kernel (—A)“°’ as K^. We start with the smooth classical solution to 
the elliptic problem 

(4.1) gu — (5div((M -|- d)'Vu) — div((u -|- dfVKc^u) = f on Bn 

with M = 0 on dBfi- 9, f ^ 0 and smooth. The sign of /, g guarantees the solution is 
nonnegative. d > 0 are constants. The nonlocal part is computed in the expected 
way by extending rt = 0 on Bf^. To find such a solution we first consider the linear 
problem 

gu — (5div((v -|- d)Vu) — div((v -|- d)VKQu) = / on Br 

for V S Cq’^ with > 0. With fixed d,5,RX >0 one can apply Schauder estimate 
theory to conclude 

with u > 0. The map T : v ^ u is then a compact map. The set {u} with w > 0 
and V S Cq’^ is a closed convex set, and hence we can apply the fixed point theorem 
(Corollary 11.2 in [11)1. to conclude there is a solution to gH). By bootstrapping 
we conclude u is smooth. 
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Now we use the existence of solutions to dnj to obtain - via recursion - solutions 
to the discretized problem 

(4.2) V^u — (5div((it + d)Vu) — div((u + d)VK^u) = / on [0,T] x Br. 

with u{0,x) = uq{x) an initially defined smooth function with compact support. 
e = T/k for some k €N. We will eventually let k —>■ oo, so that e ^ 0. 

For the next two Lemmas we will utilize the solution to 

(4.3) DTY{t) = cY{t) + h(t) 
which as in [10] is given by 


Y{t) =Y{0)Ec,{ct°‘) + a 


Jo 


s)°‘)h{t) dt, 


where E^ is the Mittag-Leffler function of order a. We will utilize in the next 
Lemmas two specific instances of (14.31) . We define Yi (t) to be the solution to (14.31) 
with Y (0) = supu(0, x),c = 0,h = 2Af. We define Y 2 {t) to be the solution to (14.31) 
with c = y2(0) = 2. and h = 0. The constant C will be chosen later. 

Lemma 4.1. Let u be a solution to (1421) . Let Yi(t) he defined as above. Then 
there exists cq depending only on T,a, ||/||l~ such that if e < eo, then 


u{ej) < Ti(ej) 

Proof. Since Yi (t) is an increasing function 

VfYiej) > A-^DfY{ej). 

Depending on T and ||/||ioo, there exists eo such that if e < eo, then 

VfYiej) > A-^DfYiej) > ^DfY{ej) = p 

We use {u{t,x) — Y{t))+ as a test function. Since {u — y)+(0) = 0 it follows from 
(12.5[) and Lemma 13.41 that 


eY,i^-Y)+V^[iu-Y)+-{u-Y).]>0. 

j<k 


We define 


B(;{u,v)= / H^{x,y)[uix)-u{y)][v{x)-viy)] dx 

JR'^ 


dy, 


where H^{x, y) — IAKq. We have the identity [4] 


B^{u,v) = / / Vu{x)'S/K^v{y) dx 

JR'^ Jr'^ 


dy. 
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Then for e small enough and j > 0 

f -y(ej))+ 

j<k 

= [ eY,iu-Y)+V^[iu-Y)+-iu-Y).+Y] 

j<k 

+ eV / {u + d)V{u-Y)+Vu 

o<k->BR 

+ e^ / {u + d)V{u-Y)+VK^u 

j<k •^Br 

> [ e^^/(M-r)++ [ (u + d)x{„>y}VMV/f^u 

J Br- ^ 


3<k 




'Br j<k ^ 


+ ^'^l^ciX{u>Y}{u + d)'^,u) 


3<k 


>xe^ / /(ej»(M(ej»-T(ej))+- 

3<k->BR 


Thus (u - Y)+ = 0. 


□ 


Lemma 4.2. Let u be a solution to g21) in [0,T] X Br. Assume that 

(4.4) 0<it(0,a;) 

If A is large there exists constants 

AIO)<5o,Co depending only on R,n,a 
€o depending only on T, a, A 
C depending only on n,a-, ||u(0,a;)||icxj, ||/||loo 

such that if pt < pLo,S < SqX < Co,e < eg, and Y 2 {t) is the solution defined earlier 
with constant C given above, then 

0 < u{ej,x) < AY 2 {ej)e~^'^^, 


for any t = ej. 

Proof. As before there exists eg depending only on T, a such that for e < eg we 
have 

VfY2{ej) > ^D^Y2{ej). 

Since u is smooth and hence continuous, u < LY 2 {ej)e~^^^ for some L > A. We 
lower L > A until it touches u for the first time. Since u = 0 on OBr this cannot 
happen on the boundary. Since u is smooth this cannot happen at a point (ej, 0). 
Also, LY 2 > 2A > 2u{{),x), so this cannot occur at the initial time. We label a 
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point of touching as {tc,rc). We compute the operator in nondivergence form and 
write Kq(u) = p and use the estimates in [5] to conclude for e small enough that 

< Vfuit,) 

= (5div((u + d)VM) + div((u + d)VK(^u) + f{tc, Vc) 

= S2[LY2{tc)e-''‘^]'^ + 6dLY2{ta)e-'''^ 

- LY2{tc)e~^^drP+ {LY2{tc)e~'^‘ + d)Ap + f{tc,rc), 


where in the equation the bar above means evaluation at Vc Then using again the 
estimates from [5, for small enough C we have a universal constant M depending 
only on n, CT such that 

|5;^|,|A^| <yi(T)M. 

Now recalling also that LY 2 {tc)e~'^'^ < Yi{T) 


-C < S{2 + d)Yi(T) + MYi(T) + 1 + 


LY2it, 


< 5{2 + d)Yi{T) + 2MYiT ( 1 + -e’'‘= ) + - . 


MYi{T) + 


f{rc)e^‘^ 

LY2{tc) 


Choosing <5, d small enough the above inequality implies 


C < 4Myi(T) +4. 

If we choose now C > 2MYi{T) +4 we obtain a contradiction. We note that C will 
only depend on n,cr, ||u(0,a;)||Loo, ||/||l~. □ 


We now give some Sobolev estimates. Because we have a right hand side we 
choose to not use ln(M) as the test function. For 0 < 7 < 1, we use {u + d)'^ — d'^ 
as a test function. The function 

FU) = -^{t + dY+^ -d^t 
7 + 1 

will satisfy the conditions in Lemma 13.41 We now assume m is a solution to (14.21) 
with assumptions as in Lemma [4.21 so that |m| < for some large M. As 

discussed in the introduction we can extend u{ej, x) = tt(0, x) for j < 0, and u will 
be a solution to (|13) on (—oo,T) X R" with right hand side 

(5div((u(0, x) + d)'Vu{0,x)) + div((u(0,a:) + d)V(—A)“'^ m(0,x)). 

for j < 0. This right hand side is not necessarily nonnegative; however, we only 
required the nonnegativity of the right hand side to guarantee that our solution is 
nonnegative. In this case we already know our solution is nonnegative. We fix a 
smooth cut-off with > M for t < —2 and 4’{t) = 0 for t > —1. We now 
take our test function as eF'{[u{t,x) — (j){t)]+). We define 

u = {u - (j))+ - {u - (j))- + (j) =: + (j). 
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We define u = u{t) for ej — 1 < t < ej. From Lemma [3.41 and the estimates in 
Section [3l there exist two constans c,C depending on a,T,A such that for e < 1, 


e^F'(u+(ej))PXej) 

3<k 


> C 




[t - s)i+“ 


"i-oo {T-tr 


dt 


-C r F'{U+{t))Dtcj,{t) dt 

J —oo 


We now consider the nonlocal spatial term. We will also use the following property: 
For an increasing function V and a constant I 


BciV{{u - l)+), u) > Bc{V{{u - l)+), {u - 0+) > 0. 
We have for the nonlocal spatial terms 



= + dy+\u) 

+ y +dr,u) 


- 7^ J ^ + d)'^+\u) 

From ProDOsition ll0.il if u'^{x) — u'^{y) > 0, then 

(uj + d)"^+\x) - (uj + d)"^+\y) > (u'^ix) - 


Then 


/ '^F'{u+{ej,x)){u+d)VKi;u> 

J<k 



Hq{x, y)\u{x)—u{y)\'^^'^ dx dy dt. 


For the local spatial term we have 


e 



"^F'{u^{<^j,x)){u + d)Vu > 7 



{u + dy\Vu\^. 
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Now combining the previous estimates with the right hand side term / we have for 
a certain constant C depending on n, a, a, A, 7 , M, T which can change line by line. 

Sj f [ (u + 


/ 0 J Br 

c'y 


7+1 


'0 

rT rt 


/ HQ{x^y)\u{x) — dx dy dt 

Jr'^ 

u{t) — u{s)]‘^ 


'BrJo Jo 
r pT 


' Br J 0 


(t - s)i+“ 

nm) 

(T - 


dt 


<C / F'{u+{t))D'^m dt 
JBr J-2 

pT p 


< 


< 


cT [ («++dr-d" 

J-2JBr 

C [ [ +dy -d? <C ( [ < C. 

J-2Jbt^ J-2 JBr 


-2 J Br 
pT 


'—2JBr J—2JBr 

The second to last inequality comes from Proposition IIP. 2( The value C is inde¬ 
pendent of C, d, R, e, d if 

(,e,6,d < 1 and R> 1. 

Then as C, d —>■ 0 we have uniform control and obtain the estimate 


Sj 


(4.5) 


'0 JBr 


(w-bd)'^lVMl' 


12+7 

llV(2-2<T)/(2+^),2+7(■s^) 


'Br 


1^11 2(0,T) — ^ 


Notice that the constant C only depends on the exponential decay of /, uq and on 
CT, a, n, T, but not on i?, 5. Letting d, C —)• 0 we obtain 

(4.6) — (5div(uVM) — div(uV(—A)”*^!!) = / on [0, T] x Br. 

We now give a compactness result. 


Lemma 4.3. Assume for any v G F, 

pT 


(4.7) 


||?;(t,a;) 


/o 


|2+7 

I14/(2-2.t)/(2+7),2+7(5^) 


®)IIw“/2.2(o,t) a C - 


Then F is totally bounded in L^’([0,T] x Br) for 1 < p < 2. 


Proof. We utilize the proof provided in for compactness in fractional Sobolev 
spaces. We will show the result for p = 2, and it will follow for p < 2 since Br is 
a bounded set. We divide T into k increments. (This k is unrelated to the number 
k for the e approximations). Let I = T/k. We define 

vi{x,t) := 7 / v{x,.s) ds. 

' Jji 
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From [9] 

(4.8) [ [ [vi{x,t) - v{x,t)] 

J Br J 0 

The above estimate is uniform for any Vj. We now utilize that [0,T] is a finite 
measure space as well as Minkowski’s inequality: the norm of the sum is less than 
or equal to the sum of the norm. 


< 


'Br 


|^c/2,2(o,T) ^ CT 


fc-1 




1 V ( 2 - 2 < t )/(2 + ',),2 + 7 ( 5 ^) 


j=o "'b 

^ 1 f /-'(J+i) \ 

•, t)||vU(2-2<^)/(2+^).2+7(B„) J 

fc -1 / 

j=0 \ 


fc -1 

= ^E 

i=o 


Jij 

■Ki+i) 


2+7 


2+7 


14/(2-2<t)/(2+7),2+7(b^)^ 


'Ij 


j{x,t) 


2+7 


W ( 2 - 2^)/(2 + 7),2 + t ( Bh ) 


It then follows from the result in [9] that for every j and A > 0 there exists 
finitely many {j3i,..., Pmj} such that for any fixed j and v G J- there exists 
I3i G {Pi,..., PMj } such that 

1 /-di+i) ^ 

v{x, t) 


IBr 


I 


'ij 


< A. 


Then combining the above estimate with (14.81) we obtain that 


/ rk-/ 3 + 4 f < / r\v-vif+[ f \vi - 

J B R J 0 J B R J 0 J B R J 0 

= [ r\v-v,\^+ij:[ b-A, 

JBrJo .nJBn 


' Br Jo 

cr + TX. 


Since I, A can be chosen arbitrarily small, J- is totally bounded. 


□ 


The following result will guarantee that V(—A) '’'u G as d —>■ 0, i? —>• oo. 

Lemma 4.4. Let u be a solution to (Irel) with right hand side f and uq both 
satisfying the exponential bound dm). Then 

f ll(-A)“'"u(t, •)II^(I-2.)/(2+7) + 2..2+^ dt<C 
Jo 

with the constant C depending only on the exponential bounds in dm. ri,^,T. 

Proof, u is extended to be zero outside of Bn. The proof is a consequence of the 
following results found in m- 

= LP{R^) n 














PARABOLIC PROBLEM 


15 


We also have the lifting property of the Riesz potential for the homogeneous Besov 
spaces 

\\{-A)-'^u\\gP+2, < C\\u\\^p^^ 

To bound u in the nonhomogeneous Besov space we recall 

||(-A)“‘^u|| 2 ,n<,/(n- 2 .,,)(Rn) < C11 U 11 i, (Rn ) . 

for any 1 < g < nj{2a). From the exponential bounds (14.41) and growth we have 
that u is uniformly in for all 1 < 9 < oo. Letting 

_ — (2 + 7)^ —^ > 1 for cr < 1/2 and n > 2, 

71 + 2(7(2 + 7) ' 

(or let (T < 1/4 for 71 = 1), we obtain by the finite length of T 

Using again the characterization of homogeneous besov spaces we obtain the result. 

□ 


Corollary 4.5. Let Uk be a sequence of solutions to (14.61) with i? —> (X) and (5—^0. 
For fixed p > 0, there exists a subsequence and limit with 

Uk ^ uq G LP{Bp) for 1 <p <2 and Uk ^ uq G 1 fU-2ct)/(2 +7),2+7^ 

Furthermore, for any compactly supported (j) 

(4.9) / [</(a^+/)2^“wo(ej,a^)+ moV(/)V(-A)"‘^mo] = eV' / /</ 

j<k j<k 

Proof. The strong and weak convergence is an immediate result of the bound (|4.5I1 
and Lemma 14.31 For 7 small enough depending on a, then 

2-2(7 


2(7 > 1. 


2 + 7 

Then from Lemma 14.41 we have that 

\/{-A)-'^Uk V(-A)-'"mo e 1F(^-2‘^)A2+7)+2^-1.2+7^ 

And in particular 

(4.10) V(-A)-'"wfc ^ V(-A)-'"mo G 

Then it is immediate from the weak and strong convergence that uq is a solution. □ 
We now show the 


Proof of Theorem 11.11 We first assume /, uq smooth and satisfying the exponential 
bounds (|4.4I) . Consider solutions Ue to (14.91) over a finite interval (0,r). As before, 
as e ^ 0 there exists a subsequence and a limit Uj —>■ Uq with the weak conver¬ 
gence as in H4.10|) and strong convergence over compact sets for 1 < p < 2 just 
as in Lemma SSI Then for fixed (p G C(j“, that uq is a solution follows from this 
convergence. The spatial piece and right hand side is straightforward to show, and 
the nonlocal time piece is taken care of as in [1]. We now consider a sequence of 
solution {uj} with {fj}, {(Tto)^} G C°° with fj —>• / and {uo)j —>■ uq in weak* L°°. 
Then again there exists a limit solution u with right hand side /. From Remark 
o and Lemma [4+] we can let T —00. □ 
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Remark 4.6. In this Section we have shown how the estimates in [5] work for 
equations of the form (12.411 . In the same way one can show that the method of 
“true (exaggerated) supersolutions” as shown in for cr < 1/2 will also work to 
prove the property of finite propagation for solutions to ()2.4[) . As the main result 
of this paper is Holder regularity of solutions we will not make this presentation 
here. 


5. Continuity: Method and Lemmas 

In this Section we outline the method used to prove Holder regularity of solutions 
to (12.41) . We follow the method used in [4] which is an adaption of the ideas originally 
used by De Giorgi. We prove a decrease in oscillation on smaller cylinders and then 
utilize the scaling property that if u is a solution to (12.4(1 , then v{t, x) = A{Bt, Cx) 
is also a solution to (12.411 if A = . Because of the degenerate nature of 

the problem the decrease in oscillation will only occur from above. Since we do not 
have a decrease in oscillation from below we will need a Lemma that says in essence 
that if the solution u is above 1/2 on most of the space time, then u is a distance 
from zero on a smaller cylinder. To prove the Lemmas in this section we will use 
energy methods, and thus we will want to use as a test function Fiu) for some F. 
If u is a solution to (12.411 . then 


g p{7(2-2cr)/(2+7),2+7 


and it is not clear that VF{u) will be a valid test function. We therefore prove the 
Lemmas for the approximate problems 


(5.1) V^u — Sdiv{D{u)\7u) — div{u\7{—A) = f on Br, 


for some large R > 0 and small <5 > 0 with u = 0 on dBji. It is actually only 
necessary to prove the energy inequalities that we will utilize with constants uni¬ 
form as i5 —> 0 and i? —>■ oo. We could also prove the Lemmas for the approximate 
problems (14.611 : however, for notational convenience and to make the proofs more 
transparent we have chosen to let e —>■ 0. Because our solution is a limit of dis¬ 
cretized solutions we then are allowed to make the formal computations involved 
with 'D^u even though u may not be regular enough for 'D^u to be defined. One 
simply proves the energy inequalities (and hence the Lemmas) for the discretized 
solutions as was done in [1]. 

Because of the one-sided nature of our problem we prove the Lemmas for solu¬ 
tions to the equation with the modified term div(ZI(u)V(—A)”*^!!), where D(u) = 
diu + c?2. We assume 0 < di,d 2 <2 and either di = 1 or fi2 > 1/2. As will be seen 
later, when d 2 > 1/2, the proofs are simpler because the problem is no longer degen¬ 
erate. We now define the exact class of solutions for which we prove the Lemmas of 
this section, u is a solution if u = 0 on dBr and for every (j) G C“((—oo, T) x Br), 
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we have 


(5.2) 


/ / / [u{t,x) — u{s,x)][(f>{t,x) — (j){s,x)]K{t,s) ds dt dx 

/ Br J —oo J —oo 

p pT p‘2t — T 

+ / / u(t,x)(j){t,x)K{t,s) ds dt dx 

Br j — CO J —CO 

— f f u{t,x)'D'^(j){t,x) dt dx 

J Br j — oo 

+ f j V(l){t,x)D(u)'S/u{t,x) dx dt 

J —CO J Br 

+ f j \7(j){t,x)D{u)V{—A)~'^udxdt 

J —CO J Br 

= f [ f{t,x)(j){t,x), 

J —CO J Br 


By Lemmas l4.3l and l4^ the Lemmas stated in this section will be true when R ^ oo 
and 5 —>■ 0. 

Before stating the Lemmas we define the following function for small 0 < r < 1/4. 
'^{x,t) := 1 + (|a:|'^ — 2)+ + — 2)+. 

We now state the Lemmas we will need. 

Lemma 5.1. Let u be a solution to (lO) with R > 4 and assume 

1 — lit < M < 'll for T < Tq 

Given /lo € (0,1/2) and tq < 1/4, there exists k > 0 depending on /ip, tq, cr,a,n 
such that if 

|{n> i/2}nr4| > (1-K)|r4| 

then u > fiQ on the smaller cylinder Li. 

We have a similar Lemma from above 

Lemma 5.2. Under the same assumptions as Lemma \5.1l given pLi C (0,1/2) and 
Tq < 1/4, there exists k > 0 depending on ^i,ro,cr, a,n such that if 

\{u>i/ 2 }nr 2 \<K\T 2 \ 


then u < 1 — fj,i on the smaller cylinder Li. 

Lemma rs.21 is not sufficient. We need the stronger 

Lemma 5.3. Under the same assumptions as Lemma \5.1l assume further for fixed 
ko 

(5.3) |{m< i/2}nr4| > (i-Ko)|r4|, 


then u <1 — pL 2 on Li for some p ,2 depending on kq. 
We will choose kq to equal the n in Lemma l5.II 
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6. PULL-UP 


In this section we provide the proof of Lemma 15.11 This Lemma is the most 
technical to prove. We first prove the Lemma in the most difficult case when 
D{u) = u + d with 0 < d < 2. Afterwards, we show how the proof is much simpler 
when D{u) = diu + d 2 with d 2 > 1/2 and 0 < di < 2. 

We will need the following technical Lemma. The proof is found in the appendix. 


Lemma 6.1. Let u, (j) be two functions such that 0 < u < (j) < 1. Let 0 < 7 < 1 be 
a constant. If |u(a;) — u{y)\ > 4|(/(a;) — <fi(y)\, then 


( 6 . 1 ) , 

Also, if 


2 /4 




,7+1 


{y) 


,7+1 


(x) 


0 < 


(/T'(y) 


< 


14 


(2;) -u^ {y)\- 


then 0 < uf,{x) - u^{y). 

If instead we assume |u(a;) — u{y)\ < 4|(/)(x) — </(j/)|, then 


( 6 . 2 ) 


,,7+1 


(y) 


,,7+1 


(x) 


< 14|(/(x) - (/(i/)|. 


(fy{y) (fy{x) 

Remark 6.2. When 0 < m < </ < 3, Lemma EH will hold with new constants by 
applying the Lemma to 


We will use a sequence of cut-off functions {</>*} which will be chosen to be 
smooth cut-off functions in space, and smooth increasing cut-off functions in time. 
We recall that for small 0 < r < 1/2, 

'^{x,t) := 1 -I- (|a:|'^ — 2)+ -|- — 2)+. 


We now recall the construction of a sequence of smooth radial cut-offs Ok from [3] 
that satisfy 

• Ok{x) < Ok-i{x) < ...< 0o{x), 

• |V6»fc|/6lfe < with m > 2 

• 9k-i — Ok > {1 — /(to)2“^ in the support of Ok, 

• Ok ^ M0XB2 as fc -)• 00, 

• the support of Ok is contained in the set where Ok-i achieves its maximum. 
We also have dp = 1 on and the support of Oq is contained in B 4 . As a cut-off 
in time we consider a sequence {^j,} satisfying 

• Ut)<^k-i{t), 

• Ckit)<c\ 

• ik->- X{t>- 2 } as fc ^ 00, 

• £,k = max^fe = 1 on the interval [—2 — 2“^,0]. 

• the support of ^k is contained in the set where ^k-i achieves its maximum. 
We now define 

(l)kix,t) := 1 - ^{x,t) + -^kit)0k{x). 

We use the convention for negative part that u = u+ — U-. We also write := 
{u — 4>k)-- We now consider the convex function 


1 

7 -I- 1 


1 

7 -I- 1 


(6.3) 


F{x) : 


(1 — xy^^ -|- a; — 
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Because of the degenerate nature of our equation we will want to utilize the test 
function 


(6.4) 


-F\u-J{(t> + d)) 


(u - ^)- 

(/>+ d 



/u + d\ 

[jTdJ 


7 



Proof of Lemma \5.1[ First Step: Obtaining an energy in time. We note that 
for 0 < X < 1, F{x) is convex,F'(x) > 0, and F"(x) > 7 . From the convexity and 
second derivative estimate we also conclude for 0 < x, ?/ < 1 


(6.5) 

( 6 . 6 ) 


F'(x)(x -y)> F{x) - F{y) + ( 7 / 2 )(x - y)"^ 
F{x) ~ x^ 


We now consider —F'{u^^/{(j)k + d))'D‘fu, and rewrite u = —u^^+cfk- To obtain 

an energy in time we first consider 


!° r p (TdlL' 

J-ooJ-00 


'^ 4 , W (s)J K{t,s) 

jt) \ uf{t)-uf{s) 

d j ^(0 ^ 


K{t,s) 


> 


rO rO 


' —00 J —00 


/ O t *0 

/ + d) 

-00 J —00 
rO rO 


(i) 


U 4 (t) u , (s) 


F 




4 >{t) + d (p{s) + d 


K{t,s) 


-F 


+ 


’ —CO —00 


{(fit)+d)^ 


+ d 

^4 (^) '^4 (' S ) 


4 >{s) + d j 

-\ 2 

K{t,s) 


K(t,s) 


/ O pO 

/ m)+d) F 

■CO J —CO 

"/“ f il 

= (l) + ( 2 ). 


+ d ())(s) -I- d 
\ +d \ (l3{s) + d 


K{t,s) 


U 4 {t)-u^{s) K{t,s) 


In the first inequality we used that (j) is increasing in t and positive for t > —4 as 
well as uf{s) = 0 for s < —4, and in the second inequality we used (16.51) . Term (2) 
is half of what we will need for the Sobolev embedding (see Lemma 12.21) . To gain 
the other half we consider term (1). For c, depending on A, a and the Lipschitz 
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constant of we have 


rO /•O 


(0fc(t) + d) 


-oo J —oo 


F 




it) 


.0 

f 

J —CO ^ 

/ —OO 


{(j)k{t) +d)F 


(l^kit) + dj 


-F 


M) 


4>k{t) + d 


(.0 nO 


+ 


l<f>k(s)-<f>k(t)]Fl 
> c / + d)F 

J —OO 


Uj,. (s) 


' —OO J —OO 

rO 


ycjikit) + d^ 

' \ 1 
^fc(0 Fd) (0 - <)“ 


(/>fc(s) + dy 

- {(t>k{s) + (i)F 

K{t,s) 


K{t,s) 

(l)k{s)+d 


K{t,s) 


rO 


> c 


C' / X{u(t)<(/ifc(i)} dt 

J —OO 


i-2-2->‘ + d)(o - ty 


C' / X{«(t)<0j.(t)} dt 

J —OO 


The last inequality coming from (16.61) . 
Now 


> 0, 


4>{t) + d ^ 

and in this proof we ignore this term which will be on the left hand side. Now for 
the term involving (pk we have 


F‘ 


J 


4>k{t) + d 


Vf^kit) > -c’^ r X{u<M- 

J —OO 


Then utilizing the embedding theorem for fractional Sobolev spaces [9] combined 

with the above inequalities we obtain 

(6.7) 

^ Uj. 


F' 


^4>k 


4>k F d 


V^u{t) > c 


o 

to 


1 

1 

1 

CM 


f f 


J — 2—2-^ J — 2 — 2-’’ 

+ 

1 


~ C / X{u<(j>k} 

J —oo 

•0 2 


_^k 


X{u<<j>k} dt. 
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After integrating in the spatial variable we have 

/ ^ \ l-ot 


rO 


{u - 4>k)^ 


{u + d)VF' 


1-00 JR" 


(pk + d 


V(-A)-'"u 


cO 


< - 


fF' 


(pk F d J 

< C"‘ / / X{u<M 


Ck 


X{u<ci>k} dt dx 


-OO . 

fO 


Second Step: Obtaining an energy in space. We now turn our attention 
to the elliptic portion of the problem. We recall from [J] the identity 

- v{y)][w{x) - ■u;(y)] 


B{v,w) = f f 
JR" JR" 


I f 

-dxdy = Cn,a / VuV(—A) '^■u. 

jR" 


\X - y|"+2-2o- 

We multiply by our test function (16.41) and integrate by parts. On the left hand 
side of the equation we have 

'u + d^'^ 


X{u<4>}^ 
X{u<(p}'y 
X{u<(p} 


(j)-\- d 
u 


-1 


(M + d)V(—A) 


(p + d 


X{u<i/>} 


7+1 


7+1 


u - 


V((m + d)/i(p + d)){u + d)V(-A)-‘"M 
dV^^ 


(p + d) 
{u + 


- 1 


((/) + d)V(-A)-‘"u 


~X{u<4>} 


{(p + d)'* 
u + dy+^ 
cP + dJ 


-{(p + d) 


- 1 


V(-A)-‘"u 

V(/)V(-A)-'"u 


7 + 1" 

:= (l) + (2). 

We now focus on (1) which will give us the energy term we need. For the term 
{-A)~'^u, we rewrite u= {u- (pk)+ - {u- (pk)- + (pk ■= + (pk- Then 

rewrite (1) = (la) + {lb) + (Ic). We focus on the term (16). We rewrite 

(16) = (16i) + (16ii) 


we 


■— X{u<(p}' 


-V 


{u + 


7 + 1 [('/' + d)'^ 


V(-A)-'^al 


+ X{«<0}:^;7PYV()>V(-A) '^ u ^. 

The term (16z) will give us the energy term in space that we will need. 


(16*) = f f 

♦/M”' «/!R^ 


{u + 

7 + 1 [('/' + d)'^ 


{x) 


X - 


{y)dxdy 


7+1 


B{X{u<4,}{u + d)^+V {(P + d)^, -u^). 
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We define the set 


Ak ■■= {|u{a:) - u{y)\ > A\(j)k{x) - (j)k{y)\}- 


It is clear that Ak contains the set 14 x 14 where we define 14 as the set on which 
6k achieves its maximum. From Lemma I6J] and Remark 16.21 we have 





(M + d)->'+^(y) 
((ji + d)T'(y) 


{u + {x) 1 (x) - (y)] (x) - {y)\^+^ 

(0 + d)'l(x) j |x-y|"+2-2- \x-y\n+2-2a 

Ak 


We now label Uk as the set where (pk achieves its maximum. Notice that Uk = 
[—2 — 2“^, 0] X 14. To utilize the fractional Sobolev embedding on 14 x 14, we also 
will need an norm of on 14. We utilize half of the integral of that we 
gained from the fractional time term: 

f r 

JKX- J-2-2-'^ (0 “ t)°‘ 

f[ ('‘iAm , 1 

- 2 JJ (0 - i)” 2 77 (0 - t)“ ■ 

Uk Uk 


The inequality comes from the fact that 0 < <1- Now from the fractional 

sobolev embedding [^, 


( 6 . 8 ) 



14x14 


\x - y|"+2-2<T 


Uk 


^ c 


n,cr,7 




(^- ^n(2+7)/(n-2+2cr) 


(n— 2 + 2 tT)/n 


This is the helpful spatial term on the left hand side that we will return to later. 

Third Step: Bounding the remaining terms. We will now show that 
everything left in our equation can be bounded by 

(6-9) [ [ X{u«i>k} dx. 

J-oc JB" 


We will denote 

A^k{x,y) := X{u{x)<<pk{x)] ~^{Ay)<4>k{v)} 

For the remainder of term (16i) we have 


in 

At 

{u + d)^+^{y) 

{u + 

Kk^^)-^^k^y)\ 

{(j) + d)^{y) 

{(j) + d)-<{x) _ 

\X - y|«+2-2o- 

r f f r ( \\Mx)-Myn 


< 

/ X{u<y>fc}- 




The last inequality is due to the Lipschitz constant of cpk when x, y are close, and 
the tail growth of cpk when x,y are far apart. 
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We now control the term (Ibii). Again, we split the region of integration over 
Ak and A%. Using Holder’s inequality (provided 2a > 7/(1 + 7) and therefore we 
must choose 7 small when cr is small) as well as the Lipschitz and sup bounds on 
(j^k we have 


( 162 ?) — Cn^(j 




- y\ 


n+2 —2(T 


dxdy. 


< r] 




- y\ 


n+2 —2(7 


dxdy. 


+ C 


< V 


[cj)k{x) - +( 2 /)]^^+^^/^^+'^^ 
k - ?/|"+2-2o- 


Xk{x,y)dxdy. 




\x - y|"+ 2 - 2 <T 

A C I \{u<4'k} 

jR” 


dxdy. 


The first term is absorbed into the left hand side and the second term is controlled 
exactly as before. 

We now consider the integration over 


<“> = // 


[Mx) - My)][u^k^x) - u^^iy)] 


\x - y\ 


n+2—2(7 


dxdy. 


< 


Ak 


[+(a;) - +(y)]2 ^ ^ 

\^_y\n+2-2a y)dxdy. 


AC / X{u<4>k} ^x. 

JR-^ 


We now turn our attention to the term (Ic). By Lemma | 6 .II we have 


|S(x{„<0}((u + d)'>'+V(+ +dy - +),+)| 

< |S(x{u<0}(u + d)^+V(+ + dy,(l)k)\ + |S(x{„<0}+),+) 


Both of the above terms are handled exactly as before by using Lemma 16.11 and 
splitting the region of integration over Ak and A|. 

The term (la) is 

( 6 . 10 ) 

(la) = S(x{„< 0 ,}u'^+V<?i'fc - (+ +d),u+J 


= 2 


X{u{x)<(t>k{x)} 


(j)k[x) + d 


(+ + d)T'(a;)_ |a; - y|”+ 2 - 2 (T 


dx dy > 0 


The factor of 2 comes form the symmetry of the kernel. We will utilize this non¬ 
negative term shortly. 

We now consider the term (2) which we recall as 


-[ [ -^-rX{u<Mii('^ + d)/i(l)k+d)y^^-l)V(j3k'^L{x-y)[u{y)-uix)]dxdy. 

jR" jR" 7 + -*- 
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In the above L = V(—A) and we have 

VL{x-y) « |x-yr(”+i- 2 <T)^ 

We again write u = + (pk- To control the term involving (pk we integrate 

over the two sets {|x — y| < 8} and {|x — y\ > 8}. We use that |</)fe(a;) — (j)k{y)\ A 
C^\x — y\ when \x — y\ <8 and \(j)k{x) — (l^kiy)] <\x — yY when \x — y\ > 8 as well 
as the bound |V(/)fc| < to obtain 


<C'" ff X{«<0,}|a;-2/1“^"”^'^^ dy 

|a:-y|<8 

+ 0 ^^ ff dxdy 

\x-y\>8 

AC' / X{u<(j>k} dx. 

jR" 


y)Y^k{y) - (l>k{x)\ dx dy 


We now use the same set decomposition with —u^^, the inequality \u^Y\ < 1 as well 
as Holder’s inequality 


7 


7+1 


X{u<M dx dy 


<CM/ dx dy 

\x-y\<8 

K,( 2 /) -U0,(a;)]2+'^ 


+ c 


\x - y\ 


n +2 — 2(7 


dx dy 


\x-v\<8 

+ ff X{u<m\^ - dx dy 

\x-y\>8 

The third term is bounded by 


C 


X{u<<j>k} dx 


provided t < 1 — 2a as well as the first term provided again that 2 (t > 7/(1 + 7 ). 
The second term can be bounded as before by splitting the region of integration 
over Ak and A'j. and absorbing the region over Ak into the left hand side. 

We now turn our attention to the last term involving . We first remark that 
the integral becomes 


X{u(x)«j>k(x)}{{{u + d)/{(j>k + - l)V^fcVL(x - y)ul^{y) dx dy. 
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We first consider the set |a; — ?/| > 8. Since < dt, 


Jj X{u(x)<4>kix)}{{{u + d)/{(l)k + - l)V(/)feVL(a; - y)u'^^{y) dx dy 


|x-y|>8 


< c 


ff X{uix)«j>u^)}\x - y\ dx dy 

\x-y\>S 

— ^ / X{u<4>k} dx. 

When \x — y\ < 8, we make the further decomposition 


(l>k{x) 


-y\<v 


to absorb the integral by the nonnegative quantity (16.101) . In the complement when 

NMx)\, 


(l>k{x) 


-y\>v 


we use (jik > \V(j)k\/(l>k and integrate in y 


VL{x - y)ul^{y)dy 


y,n—l 


< 




^ < max{C, 


The remainder of the terms are bounded by |V0fe| < By multiplying by 

the term X{u<ij>} and integrating, we end up in the worst case with 

/^k f — l/77i+(2cr— 1)/771 7 ^ /^k f J 

C I X{u<(pk}4^k dx ^ C j X{u<^k} dx^ 

J'MX 

since m > 2. The last term to consider is the local spatial term. We use Cauchy- 
Schwarz 


d / X{u<<i>k} 

Jbr 


V 


u + d 


7 


= 5 7 


\4'k + d 
u + d y |Vup 


- 1 


(u + d)VM 


Bjt \4'k + dy (pk + d 


-cs 


7+2 


X{u<^k} 
7 


IBr \(t>k + dy 

- <' ■ {BB) 
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Retaining the energy from (j6.8l) on the left hand side and moving everything else 
to the right hand side which is bounded by (16.9p . our energy inequality becomes 


c { 

JVk V“'-2-2-'= 


— ',n( 2 + 7 )/(n—2+2 ct) 


(n—2+2cr) / n 


<C^ f f X{u<M<C'^ 

JM.'^ J —oo J J 

Uk-i 

Fourth Step: The nonlinear recursion relation. We now (as in m) use 

Holder’s inequality twice with the relations 

P . 1-/3 1-/3 


Pi P2 P P3 Pi 


for a function v to obtain 

yP < f [ 


< 


< 


V- 


pP/Pl / n \ p(1-/3)/P2 

-Pi 1 ( / .,P2 


Pa/pi 


I V- 
0pIP3 


We now choose 


Pi =2, P2 = 


n (2 + 7 ) 
n — (2 — 2a) ’ 


■ X P4/P2\ 

n 


P3 = -. -, P4 = 2 + 7, 

1 — a 


so that if r = 2 — 2(7 


r + an (2 + 7)/2 r 

P = 2 — 7 ;-^- and p = 


(1 — a)r + an r + an(2 + 7)/2 

We now use Holder’s inequality one more time to obtain 


hjp 

yP 1 < 


nPl 


P3lpi\^^/P^ 


1 

< — 

OJ 


Pa/pi 


Pbu>IP3 


„P2 


OJ — 1 


(jJ 


P 4 /P 2 


{l-0)b/pi 


(1-I3)buj 

• \ P4/P2\ P4(‘^-l) 

j 


We choose 


(2 + 7/3) ^ ^ 2(2 + 7 ) . 

and 0 = < p 


(2 + 7)/3 


2 + 7/3 


UJ = 
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where we used Minkowski’s inequality in the last inequality. Substituting for 
V in (|6.12p and utilizing (16.8p we obtain 


(6.13) 


/ N 


//(»-«' 

II X{u<M 

\Uk / 

' Uk-i 


We first recall that > (pk + — /ro)2 ^. We now utilize Tchebychev’s 

inequality 

Jf X{u<M < (2V(1 - JJ{n- 

Uk Uk-i 

Combining the above inequality with (16.131) we conclude 


If we define 



( \ 

JJ{u- PkY- < 

Iliu-Pk-iY 

Uk 



Mk ■■= JJ{u - pk)-, 

Uk 


Then 

Mk < 


Since p > b, ii Mq is sufficiently small - depending on C and p/h - we obtain that 
Uk —>■ 0, and hence u > po. 

□ 


We now prove Lemma l5.II in the case when D{x) = {dix + ^ 2 ) with di < 2 and 
^2 > 1/2. This is actually much simpler because we can utilize the test function 
— {u — (pk)- as when dealing with a linear equation. 

Proof. We choose as the test function —{u — pk)-- Notice that 
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The fact that d 2 > 1/2 gives a nondegenerate linear term which we utilize. From 
the computations in [1] we then have 

/ r 


>C 

Jl 

_Ck 


•oo ^ ^ 

/ / X{u<0fc} 

J — oo 


and even more importantly 
^0 


(6.14) 


/ / Vir^^V(-A)--u 

J —OO c / r ^ 

[ [u-^S^)-u-^{y)Y 

-oo Jr" |a: - y|’"+2-2'^ 

/ / X{«<0,} dx dt. 

J-oo JR" 


> C 


dx dy dt 


-C* 


Notice that in the above inequality we have the power | • p rather than | • p+T". 
We now show how to bound the terms invovling diu. 

~ %k + = V[(u"jV2 - u-^(j)k] + u-ycjik = (1) + (2) 

Then multiplying (1) by V(—A)~°'u and integrating over R” we have 

Now 

[(“0JV2 - - [(W0JV2 - u^^<pk]iy) 

= ^ (y) - ( 2 ;)] [</fe (x) + (t)k (y) - u-^ (a;) - u-^ (y)] 

- \ iy) + (a;)] [</fc (x) - <l)k (y)] 

= (la) + (16) 

We write u = — Uipf.. We break up our set into the two regions 

Fk ■■= {{x,y) : \u~{x) -u~{y)\ > 2|(/(x) -<(>(y)|} 

We notice that on the set F we have that 

(l){x) + cj){y) - u^{x) -u^{y) > \u^{x) -u^{y)\/2. 

Then integrating over F we have for the term (la) with right term — 

/2 - u-Jk]ix) - [(u^jV2-u“ (/fe](i/))(u^-(^) -M“ (y)) 


Fk 


' - y\ 


n+2—2(T 


> 


Fk 


\x-y\ 


n+2 —2f7 — 


> 0 . 
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This is the nonnegative energy piece which we actually do not need having obtained 
a better piece in (j6.14ll . All of the remaining terms in (1) can be bounded by break¬ 
ing up the region of integration over Fk,F^. Over Fk we use Holder’s inequality 
with p = 2 rather than with p = 2 + ^ and absorb the small pieces by the term 
in (16.141) . We use the same methods as before to bound the integration over F^. 
Bounding the term (2) is done as before with slightly easier computations. 

The local spatial term is bounded in the usual manner. □ 


7. Pull-down 

In this section we prove Lemma [521 We will need the following estimate that is 
analogous to Lemma [6.II 


Lemma 7.1. Let u,(f> be two functions such that 1/2 < /> < u < 1. Let 0 < 7 < 1 
be a constant. If |u(a;) — u{y)\ > 4|(/)(x) — />(2/)|, then 


(7.1) 

Ci\u^{x) - u+{y)\^+'^ < 

mT'+1(?/) 

u'^'^^{x) 

(tA{y) 

(j)'y{x) 

Also, if 





n < 

u'^~^^{x) 



<l>^{y) 

^T'(x) 




then 0 < {u — 4>) - (x) — {u — 4>) - {y). 

If instead we assume |u(a;) — m(?/)| < 4|(/)(x) — (f{y)\, then 


(7.2) 




<l)'f{x) 


< lA\(j){x) 




The proof is similar to the proof of Lemma 16.11 In this case since u> (f one uses 
the bound above on u and the fact that (f> is bounded by below. 


Proof of Lem,ma \5.Si The proof is nearly identical. We mention the differences. We 
only consider D{u) = u since the modifications for handling D{u) = diu + d 2 have 
already been shown in the proof of Lemma l5.ll We consider a similar test function 


F{x) 


—-(1 -I- xy^^ - X + 

j + 1 


1 

7 -I-1 


We then utilize 

(7.3) (^1+ 


This time we consider the same test functions Ok{x) in the space variable, but this 
time we multiply only by a single cut-off in time foit)- We define our ifk as 

4>k ■= t) - fo{t)0kix/2). 

To obtain the same estimate in time we only need to recognize that fik is now 
decreasing in time and bounded by below by 1/2. 
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F' 






= mF' 

= mF 


„ ( u+{t)-u+{s) 


Ht) 




(7.4) 





1 , m J 

_ ^(t) 

fi(s) _ 




4 >(t) <j){s) 


> - F{u+{t)/(l){t)) - F{u+{s)/(l){s)) 


- 

2 L 

-Cu^is)it-s) 

The negative constant comes from the fact that cj)~^ is Lipschitz. Then everything 
proceeds as before. Since our cut-off is bounded by below our norm in time 
occurs over all of (—oo,0). We obtain as before 


{u - (/>)+ 


' —OO 

rT 


uVF' V(-A)- 


< 


X{ii> 0 } dt dx 


The spatial portion of the problem is handled exactly as before. 


□ 


8. Decrease in Oscillation 

We define 

F(t, x) := i sup{-l, inf{0, \x\^ - 9}} -I- ^ sup{-l, inf{0, \tf - 9}}. 

We point out that F is Lipschitz, compactly supported in [—3,0] x B 3 and equal 
to —1/2 in [—2,0] X i? 2 . We also define for 0 < A < 1/4, 

■■= (d^l - A-i/")" - 1)+ + ((|t| - - 1)+ for |t|, |x| > 

and zero otherwise. The value of v will be determined later. Finally, we define for 
iG{0,l,2,3,4} 

</i = 1 + tpxs + X^F. 

Then 1/2 < (/o < ■ • ■ < 1/4 < 1 in r 4 . 

Lemma 8.1. Let k be the constant defined in Lemma \5.‘A Let u he a solution 
to ([5.2|) . There exists a small constant p > 0 depending only on n,a,a and Aq 
depending only on n, a, a, p, 6 such that for any solution u defined in (a, 0) x R” 
with a < —4 and 

—ijjx 3 < u{t, x) < 1 -I- ipx^ i'n (a, 0) x R" 
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with A < Ao, and / < A^, then if 

|{w < (jjo} n {Bi X (-4, -2))| > p, 


then 

|{u>(?i 4 }n(R” X (-2,0))| <K. 

Proof. We will show the computations for D{u) = u. The general situation is 
handled as before as in Lemma oi 

First Step: Revisiting the energy inequality. We return again to the 
energy inequality. This time, however, we will make use of the nonnegative terms. 
We seek to obtain a bounde on the right hand side of the form 

We now consider the test function as in (El, but with cut-off (fi- It u > (fi , 
then 1/2 < (pi < u < 1, and so 

yl — (P 

F'{u/pi) = — < 27 U+ < 27 A^\ 

To take care of the piece in time we first note that pi is Lipschitz in time for 
t € [0,4] with Lipschitz constant 2A. Then as before 


J F'{u/Pi)V’/u dt = J F\u/Pi)V/{u+^ - +P) = (Tl) + (T 2 ) + (TS). 

For (Tl) we return to the inequality (17.41) and utilize the Lipschitz nature of pf^, 
to obtain 


F'{u/pi)V/u+^ > c / {u- pi)/: 


' —00 
r*0 rt 


1-a 




> C 


— 00 J —00 
nO 


Mt)upFfuJ^{t)/p{t)){P^ it)-pi (s)) dt 

. 1-a 

{u - pi)/r 


' —00 

-pLp'- 

The nonnegative piece (T2) will be utilized in the second step of this proof. For 
(T3) we note that since pi is decreasing, we have 

0 < -V/Pi < -K-’^DfPi = -T>“i/>a3 - DfXF. 

Clearly, —DfyF < CA® for t < 0 from the Lipschitz nature of F. For —4 < t < 0 
we have 

/•A-'/" I I 




ll-l-Q 


We therefore pick v small enough that {a — v)/v > 2. 


^0 ^0 
/ F’{ulPi)V/Pi <C Xu+< CXP 


— 00 
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Our energy inequality becomes 


c / 


+ c 




Ki 


(t-s) 


l+a 


+ “Spatial Terms” <CX^+ fuj^. 


Since / < A^, everything is bounded on the right hand side by CA^. 

We now turn our attention to the elliptic portion. We consider the terms 
(la), (16i), (16ii), (Ic), (2) as the analogous terms for those defined in the proof 
of Lemma O As before we obtain a nonnegative energy from the term (16i). Ev¬ 
erything else we will absorb into this energy or bound by C ^The term 
from (Ibi) over AJ is bounded as follows: 


[(t>i{x) - ((.i(y)][u+(x) - u+(y)] 


n+2—2cr 


< c 


I \x-y\ 

,|n+2 —2(7 


dx dy 


dx dy 


I A- \x-y\‘' 

We have the following inequality from the computations given in |3]: 

( 8 . 1 ) ff - ^ 

JJ \x-y\'^+‘‘ 

with (2 — 2cr — 2v)lv > 2. In particular, (18.11) will hold for 7 = 0. For the term 
( 16 m ) we break up the region of integration into Ai and Af. On Ai we use Holder’s 
inequality as before 

[hjx) - (^1 (y)] (a;) - ^^^^(y)] 

/ All 


< C 


\x - y|"+2-2<T 

[1^1 (a:) - 

\X - y|"+2-2o- 

[u^,{x) - u^,{y)]‘^ 


dx dy 


dx dy 


Ai 


< -f y 


, , „ „ dx dy 
X — y|”“'"^ 

[u,j>r{x) -U4,^{y)Y 

Ai |a:-y|"+2-2- 


dx dy 


The last term is absorbed into the left hand side. The other term is bounded again 
from (ED. The term (Ic) is bounded in exactly the same way. (la) is nonnegative 
and will be utilized later. We now turn our attention to the term (2). We rewrite 
u = + cf>i. The term involving with |a: — y| < y is absorbed by the 

nonnegative term (la) on the left hand side. We now utilize the inequalities: 

• |VL(a: — y)| ss l/|a: — 

• jv^il < C for all X 

• |V0i| < CA in the support of 

• “ 1] ^ 4u+ < 4A. 
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Then 


/R" 7 ■ 


< C\^ 


-^X{u< 4 >}{{u/- l)V^iVL(x - y)[u{y) - M(a;)] dx dy 
X{«< 0 }VL(a: - y)[u{y) - u{x)] dx dy 


/R" ^R" 

These terms are all bounded as before. Notice that we have on the outside. 
Then all nonloacl terms on the right hand side are bounded by The 

local term div{D{u)\/u) is handled in the usual manner by use of Cauchy-Schwarz. 

Second Step: Using the “good” spatial piece. We now utilize the two 
nonnegative pieces. From Proposition IIP, ll we have 

7 + 1 


-V 


- 


>4(<) 


+ P+T 


-I + 


Then we conclude that 


( 8 . 2 ) 


pO /•o 


/-4 o 
rO 




1 + Q! 


t — s 


{u^,V+\x)u^^{y) 


/-4 . 


X 


-y\ 


n+2—2(T 




Since we used ipAa, replacing (j)i with we have the same inequality but with the 
bound CA 3 ( 2 + 7 )/(i+ 7 ). 

We now show how the inequality (18.21) and its analogue for (j)^ are enough to 
prove the remainder of the Lemma as in [1] . We note that for the proof as written 
in [1] to work we need 3(2 + 7 )/(l + 7 ) > 5 which is achieved for 7 small enough. 
We first utilize 

f f ^ (2+7)/(1 + 7) 

VV |:r-yr+2-2. - 

From our hypothesis 

|{w < 00} n ((-4, -2) X Bi)\ > p. 



Then the set of times E S (—4, —2) for which |{u(t, •) < 0o}ni?i| > p/4 has atleast 
measure p/{2\Bi\). And so 

C'A(2+t')/(i+'>') > cp / [ {u+J^+'-^dxdt 
JT. Js." ^ 


Now ({m — 02 > 0} n (S X B 2 )) C ({u — 01 > A/2} n (S X B 2 )), and so from 
Tchebychev’s inequality 

|{m - 02 > 0} n (e X ^2)1 < 

P 

The exponent on A is positive for 7 small enough. We write this as 

|{u < 02} n (E X ^2)1 > |E X S 2 I - > p/2 - 

P P 

This will be positive for A small enough depending on n, cr, a, 7 , p. 
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The proof then proceeds just as in [1] where we then utilize 3(2 + 7 )/(l + 7 ) > 5 
as well as the analogue of (| 8 . 2 I 1 for ^ 3 . □ 


This next Lemma will imply Lemma 15.31 For this next lemma we define 


lA.,. = ((|x| - - 1 )+ + m - - 1 )+ 


Lemma 8.2. Given p > 0 there exist r > 0 and pi such that for any solution to 
(EH) in R” X (a, 0) with a < —4 and \ f\ < satisfying 


If 


then 


< U < 1 + V'r,A, 


|{u < 0o} n {Bi X (-4, -2))| > p, 


sup U < 1 — Pi- 
Bix(-1,0) 


Proof We consider the rescaled function w{t,x) = (u — (1 — A'^))/A"‘. We fix r 
small enough such that 

< (Ixy/^ - 1 )+ and < (|tr/4 - 1 )+ 

Then w satisfies equation (I5.2|) with D 2 {x) = Di{X‘^x + (1 — A"^)) where Di is the 
coefficient for the equation u satisfies. From our hypothesis and Lemma l8.II 

|{u; > 1/2} n (S 2 X [- 2 , 0 ])| = |{u > M n (B 2 X [- 2 , 0 ])| < p. 

Then from Lemma [521 we conclude that w < 1 — pi on (—1,0) x Bi, and so 
M < 1 — A^pi = 1 — P2- Q 


9. Proof of Regularity 

With Lemmas isnioi and 15.31 we are ready to finish the proof of Theorem 1 1.41 
We first mention that solutions of (12.41) satisfy the following scaling property: If u 
is a solution on (a, 0) x K", then v{t, x) = Au{Bt, Cx) is a solution on {a/B, 0) x R" 
ii A = B°‘C^'^~^. The method of proof is given in s which we now briefly outline. 
We take any point p = {x, t) G R" x (a, T) and prove that u is Holder continuous 
around p. The Holder continuity exponent will depend only on a, a, n. The constant 
will depend on the norm of u, f and on the norm of u(a, x). By translation 
we assume that p = (0,0). By scaling we assume that 0 < u{t,x) < 'i>(t,x) and 
I/I < A3 for A as defined in - W 1- W 1 

We now take a positive constant M < 1/4 such that for 0 < if < M 

Aia:). 

M will depend only on A, p 2 and r > 0. During the iteration we have the following 
alternative. 

Alternative 1. Suppose that we can apply Lemma 15.31 repeatedly. We then 
consider the rescaled functions 

1 / 

—^Uj(Mit,Ma:), Mi = I 


Uj + l{t,x) 


1 - P2/4 
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Notice that Mi < M. All the Uj satisfy the same equation. If we can apply Lemma 
I5.3l at every step, then Uj < 1—^2 on the cylinder Li. This implies Holder regularity 
around p and also implies u{p) = 0 . 

Alternative 2. If at some point the assumption (lOI) fails, then we are in the 
situation of Lemma 15.11 and 

0 < ^0 < Uj{t, a;) < 1. 

Scaling the above situation our equations will have D{u) = diu + d 2 with ^2 > 0. 
We may then repeat the procedure since Lemma 15.11 and 15.31 apply also in this 
situation. 


10. APPENDIX 


Proof of Lemma Id. 11 Since throughout the paper we only require 7 small when a 
is small, we will prove the Lemma for 7 = 1 /fc for fc € N. Now we assume without 
loss of generality that 


,7+1 


{y) 


, 7+1 


{x) 


> 0 . 


4n{y) 4n{x) 

We first assume that |it(a;) — u{y)\ < 4|(^(a;) — (/>(y)|. We need to bound 


,7+1 


{y) 


,7+1 


{x) 


We first notice that term above in (110.11) will be larger if we assume that u{y) > u{x) 
and (j){x) > 4>{y) without changing \u{x) — u{y)\ and \4>{x) — 4>{y)\- Furthermore, 
the term in ( 110 . 11 ) will still be greater if u{y) = (j){y) and not changing \u{y) — u{x)\. 
We are then looking for the bound 

“( 2 ^)- A-yf \ - - ^iy)\ = C 2 ( 0 (a;) - u{y)). 

(tP{x) 

Thus, for a constant I we need the bound 

( 10 . 2 ) 


I Ap <c,. 


{I + 

Recalling that we are assuming 4:\(f){x) — (j){y)\ > \u{x) — u{y)\ (or v < 4/r) the above 
term is maximized when v is largest or when v = 4/r. Now 

(Z-4/r)i+T' 


I - 


{I + /i)^ 

i[{i+p)^ -{i-4pr] 


{I+p)^ 
■= Li + L2 


pff-M^ 

(2 + /i)-" 


It is clear that 


L 2 A Ap. 

To control Li we first consider when I — Ap < 1/2. Then I < 8p and it is clearly 
true that 

Li < 8p. 

Now when I — Ap > 1/2, from the concavity of x'^ we have 

I {I - Ap)"< 


Li A 


I — Ap (Z + pY 


-5p < lOp. 
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Then 

(10.3) I - < 14/x, 

and (lOI) is proven with constant C 2 = 14. 

We now assume 4|(^(a;) — (j){y)\ < \u{x) — u{y)\, and the left hand side of (110.211 
is maximized again when Afj, = v and so we have 

(Z-u)i+T' ^ 14^^ 

(/ + u/4)T' “ 

which is just (110.31) rewritten with the substitution /i = r;/4. Then 


(10.4) 


,7+1 


(.y) 


.,7+1 


(x) 


14 


14, 


(/>^iy) 


- u{y)\ < y (x) - y (y)|. 


and the right hand side of dnu) is shown. 
Now 


,7 + 1 


(y) 


,7 + 1 


(x) 


7+l(y)-„7+l 


{x) 






^+1, {x) - (jp {y) 

'X ' [x ]— , , , , , , := Ml + M2 

^ ’ cjx^{x)ctx^{y) 


We suppose y = Ijk. By factoring we have 

^ \u{y) - u{x)\ Ejto n("-^')/fe(y)M/fe(x) 

> u{x) _ 4|(/)(a:) - (j){y)\ 

(tp{y) X)j=o u^^~'^~^'^/^(,y)u^^^{x) 

1 ^ I ^ _ \(l){x) - (i){y)\ 

^ u{x) _ \(t>{x) - (/)(y)| 

“ 4P{y) X)j=o u^^-^-iy'^(]j)ui/’^{x) 

Thus M 2 < Ml!A. Thus, Mi is the dominant term. We then have from the 
convexity of 


Ml + M 2 > Mi/2 > 


> 


> 


> 


1 u^^^{y) — u^^^{x) 

2 

1 {u{y) -u{x)Y+'^ 

2 cjx^iy) 

, , ,1+^1 

^ ^ 2 xjx^iy) 

5^ ' > xjx^[y) 


□ 


Proposition 10.1. Let F be a function satisfying F,F" > 0 for x > 0. Assume 
also F{0) = 0. If y > X > 0, then 

F{y) - F{x) > F{y - x). 
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Proof. For fixed h > 0, 

d F{x + h) — F{x) F'{x + h) — F'{x) ^ ^ 
dx h h ~ 

Then for x,h >0 

F{x + h)- F{x) ^ F(0 + h)- F{0) _ F{h) 
h ~ h h 

Let h = y — X, and multiply both sides of the equation hy y — x. 

Proposition 10.2. Let 0 < 7 < 1. Let x,d> 0. Then 

{x + df -d? < 2^x'< 

Proof. First assume x < d. From the concavity of x'*' we have 

{x + dy — dP < ydP~^x = ydP~^x^~^x^ < yx^. 
If on the other hand x > d, then 

{x + dy -d'< < {x + dy < {2xy. 


□ 


□ 
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